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ABSTRACT 

Aims. Deflagration models poorly explain the observed diversity of SNIa. Current multidimensional simulations of SNIa predict a 
significant amount of, so far unobserved, carbon and oxygen moving at low velocities. It has been proposed that these drawbacks can 
be resolved if there is a sudden jump to a detonation (delayed detonation), but these kinds of models have been explored mainly in 
one dimension. Here we present new three-dimensional delayed detonation models in which the deflagraton-to-detonation transition 
(DDT) takes place in conditions like those favored by one-dimensional models. 

Methods. We have used a smoothed-particle-hydrodynamics code adapted to follow all the dynamical phases of the explosion, with 
algorithms devised to handle subsonic as well as supersonic combustion fronts. The starting point was a centrally ignited C-O white 
dwarf of 1.38Mo. When the average density on the flame surface reached ~ 2 - 3 x 10' g cm"^ a detonation was launched. 
Results. The detonation wave processed more than 0.3 Mq of carbon and oxygen, emptying the central regions of the ejecta of 
unburned fuel and raising its kinetic energy close to the fiducial 10^' ergs expected from a healthy Type la supernova. The final 
amount of '^Ni synthesized also was in the correct range. However, the mass of carbon and oxygen ejected is still too high. 
Conclusions. The three-dimensional delayed detonation models explored here show an improvement over pure deflagration models, 
but they still fail to coincide with basic observational constraints. However, there are many aspects of the model that are still poorly 
known (geometry of flame ignition, mechanism of DDT, properties of detonation waves traversing a mixture of fuel and ashes). 
Therefore, it will be worth pursuing its exploration to see if a good SNIa model based on the three-dimensional delayed detonation 
scenario can be obtained. 
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1. Introduction 

The knowledge of the physical mechanism by which a white 
dwarf is disrupted by a thermonuclear explosion is relevant to 
many topics of modern astrophysics. A satisfactory model of 
the explosion becomes crucial to better understand Type la su- 
pernovae, which in turn have profound implications in cosmol- 
ogy and in studies of the dynamics of the interstellar medium 
and the chemical evolution of galaxies. Nevertheless the de- 
tails of the explosion mechanism are not well known. Nowadays 
the best models involve a white dwarf with a mass near the 
Chandrasekhar-mass limit in which a nuclear fuel ignites in one 
or many sparks around the center and thermonuclear combus- 
tion propagates rapidly through the rest of the star. In these 
models the starting point of hydrodynamical calculations is the 
thermal runaway caused by the screened '^C-h'^C nuclear re- 
action in degenerate conditions. The shape of the initially in- 
cinerated region probably does not keep the spherical symme- 
try that characterizes previous, hydrostatic, stages due to the 
stochastic nature of ignition and to the role played by hy- 
drodynamic instabilities. As a consequence, any comprehen- 
sive numerical simulation of the explosion has to be performed 
in three dimensions. Many exploratory studies carried out in 
the past using one-dimensional hydrodynamics have shown 
that the key point to achieve a successful explosion relies in 
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the proper estimation of the velocity of the combustion fr ont 
("Nie mever & Wooslev 1997; Hillebrandt & Niemeyer '2000). It 
has been shown (Timmes & Woosley 1992) that the width of the 
flame is about Ifim at p ^ 10^ g cm"^, implying that at scales of 
a km the flame can be regarded as a thermal discontinuity sepa- 
rating ashes from fresh nuclear fuel. Although the flame velocity 
is large, close to 100 km s"', it is quite subsonic: ^ O.Olvsound, 
where Vjound is the sound speed. These velocities are far below 
the nearly supersonic velocities which, according to many one- 
dimensional simulations, are needed to blow away the star. It is 
today widely accepted that the hydrodynamic instabilities inher- 
ent to the explosion are the drivers of combustion acceleration, 
thanks to the huge increase in the effective area of heat exchange 
between fuel and ashes. 

Up to now the majority of three-dimensional studies of Type 
la s upernovae have focused on the p ure deflagration paradigm 
(Ga rcia-Senz et al.l 19981: Reinecke e t al. 2002; Gamezo et al] 
2003; 'Garcfa-Senz & Bravo"2005'; 'Schmi dt & Niemeveil 120061; 
Schmidt et al. 2006; Ropke et al. 2006b a). In these models the 
flame front moves fast but always subsonically, leaving time for 
the white dwarf to expand before the deflagration reaches the 
outermost regions of the star. These models have been success- 
ful in explaining many observational facts using a minimum set 
of free parameters. Unfortunately, all the three-dimensional ver- 
sions of the pure deflagration paradigm present several weak- 
nesses: 
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- The presence of a few tenths of unburnt, and unseen, car- 
bon moving at low velocity is probably the worst flaw 
([Kozma et al. 2005). 

- Although the amount of Fe-group elements synthesized in 
the explosion is within the expected range, the actual mass 
of ^''Ni ejected might turn out to be too low when neutroniza- 
tion due to electron captures is properly taken into account. 

- The final kinetic energy never reaches the canonical value of 
1 foe (1 foe = 10^' ergs). 

- The synthesis of i ntermediate-mass elements is scarce (see 
iGarcfa-Senz et al.l l2007). 

- The ejecta are not chemically stratified, as demanded by ob- 
servations of supern ovae and remnants (iBadenes et al.ll20d6t 
iGerardv eTaHbOOTh . 

- Large clumps of radioactive ^^Ni and ^^Co are present at the 
photosphere at the time of maximum brightness. 

There have been several attempts to bypass the weak- 
nesses of deflagrations models, preserving at the same time 
many of their advantages, by postulating a sudden transi- 
tion from a subsonic de flagration to a supersonic detona- 
tion dlvanova et al ] 11974 . i.e. the multidimensional version 
of the delayed deton at ion m o del: Khokhiovl (Il 99l|) in ID: 
[Arnett & Livnd (Il994l): lUviiS (Il999h; iGplombek & Niemever 
(|2005ft in 2D: Garci'a-Senz & Bravo' 7 20031) : iGamezo et al. 
( |2004 |2005): Ropke & Niemeyer (2007) in 3D. In particular, 
the recent work of iGamezo et al.l (l2005h gave a satisfactory ther- 
monuclear explosion solving many of the problems of pure de- 
flagration models. The density at which the jump to the detona- 
tion takes place had to be estimated because the physics behind 
the deflagra tion-to-detonat ion (DDT) transition is not well un- 
derstood. In lGamezo et al.l (12005) the adopted transition density 
was rather high, pt > 2.5 x 10** g cm"^ an d the detonation was 
induced near the center of the white dwarf. iRopke & Niemeverl 
(1200 7) explored the range of explosions that would result if the 
DDT took place in the outer shells, at r ~ 1.7 - 1.8 x 10** cm, 
after a huge expansion drived by subsonic combustion initiated 
in a few hot bubbles. They used a physically motivated criterion 
to switch on a DDT, namely that the combustion entered the dis- 
tributed regime, in which mixing of burnt and fresh fuel by tur- 
bulence fluctuations at the scale of the flame width might trigger 
a detonation. Although this criterion is as fragile as any other 
currently used in multidimensional simulations of DDT, it led 
to consistent results with respect to the final kinetic energy and 
to the mass of Fe-group and intermediate-mass elements synthe- 
sized. It also allowed almost complete depletion of carbon in the 
center 

Given the scarcity of three-dimensional simulations of the 
delayed detonation scenario and the complexity of the problem, 
it is interesting to perform additional calculations to explore fur- 
ther the range of initial conditions of the DDT. In this work we 
present new three-dimensional simulations of delayed detona- 
tions in white dwarfs in which the DDT takes place in conditions 
close to those usually assumed in one-dimensional calculations. 
In the present models the DDT was triggered when the mean 
density at the flame surface was p~2-3xl0^g cm"-'. We have 
also explored the sensitivity of the results to the initial geometri- 
cal configuration of the flame by computing a delayed detonation 
model in which the deflagrative phase started at the center of the 
white dwarf, and another model in which i t started in 30 hot bub- 
bles, a s in the deflagration model B30U in lGarcia-Senz & Bravol 
(I2005h . 

The plan of the paper is as follows. In Sect. 2 we briefly 
describe the hydrocode and the physics included in the three- 



dimensional simulations. Next, we discuss the hydrodynamical 
evolution during both the deflagration and the detonation stages, 
as well as the final nucleosynthesis. Some discussion of our re- 
sults and a comparison to previous works as well as the conclu- 
sions are provided in Sect. 4. 



2. Methods 

The simulations were carried out using a smoothed-particle hy- 
drodynamics (SP H) code suited to ha ndle both subsonic defla- 
gratio n waves (.Garcia-Senz et al.llT998.) and supersonic detona- 
tions (iGarcia-Senz et al.) Il999h . The calculations were carried 
out using 250000 particles, which translated into a maximum 
spatial resolution of 22 knfl (about 1 % of the white dwarf radius) 
during the deflagration phase and, later on, of 15 km (~ 0.2% of 
the stellar radius) during the detonation phase. The equatio n of 
state, the same as described in lGarcia-Senz & Bravol (l2005h . in- 
cludes all the terms relevant for the explosive phases of Type la 
supernovae. The properties of matter in nuclear statistical equi- 
librium (NSE) were interpolated from accurate tables giving the 
nuclear binding energy, electron capture and neutrino generation 
rates as a function of the temperature, density and electron mole 
number. This represented a negligible overhead of CPU time and 
allowed us to obtain realistic estimates of the nucleosynthesis 
during the SNIa explosion. As usual in this kind of simulation, 
the detailed nucleosynthesis was calculated by postprocessing 
the output of the hydrodynamical model. 

During the deflagration phase the flame was propagated by 
a diffiision-reaction scheme devised to match a prescribed flame 
velocity law. Although the SPH algorithms are flexible enough 
to incorporate a subgrid-scale model able to capture the eff'ects of 
turbulent flows into the effective flame velocity (Garcfa-Senz et 
al. 2008, in preparation), for the present calculations we choose 
to stay as close to possible to the usual assumptions of suc- 
cesful one-dimensional delayed detonation models. Thus, we 
implemented a constant flame velocity in the range Vflame = 
100 - 200 km s ' , independ ent of density, durin g t he defl agra- 
tion phase. As described in iGarcia-Senz & Bravol (l2005l) . this 
prescription has proven adequate to describe the gross features 
of subsonic burning propagation. Note that the true laminar ve- 
locity in the center of the white dwarf is close to vgame - 100 km 
s ' . On the other hand the time-averaged effective deflagra- 
tion velocity which stems from that scheme is similar to the 
one assumed in successful one-dimensional delayed detonation 
models during the deflagratio n phase , ^ 0.03vsound - 0.05vsound 
('Dommguez & Hoflich"2000). 

During the detonation phase the above described algorithm 
for flame propagation was switched off' and the combustion was 
followed by solving the nuclear kinetic equations together with 
the hydrodynamic and energy equations. The nuclear part con- 
sisted of a small nuclear n etwork of 9 nuclei from helium to 
nickel (iTimmes et al.ll20()0h with prescriptions aimed to handle 
both quasi-statistical equilibrium and nuclear statistical equilib- 
rium. Whenever the temperature became higher than 5.5 x 10^ K 
the material was assumed to be incinerated to NSE if the density 
was higher than 10^ g cm"^. Once the NSE regime was attained 
the abundances were not computed by solving the nuclear kinetic 
equations, instead the chemical composition was obtained from 
the NSE tables. In our calculations, incinerated particles were 
assumed to leave NSE when either their temperature or their 



' Here we give the smoothing length of the SPH code as a proxy for 
the spatial resolution. Note that the smoothing length is both spatially 
and temporally variable 
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density were lower than 2 x 10^ K or lO*" g cm"^, respectively. 
During this phase the detonation front was captured, that is the 
particles within the detonation front were identified with the aid 
of an algorithm that looked for simultaneous strong gradients of 
pressure and temperature together with sonic flows. Particles lo- 
cated at a detonation front were assigned ellipsoidal kernels to 
perform the SPH interpolations in such a way that the smallest 
axis of its ellipsoid was aligned with the normal to the detonation 
front. In regions where there were no shocks the kernel reduced 
to the standard spherically symmetric kernel. The ability of this 
de vice to represent detonati on waves was demonstrated in Fig. 1 
of iGarcia-Senz et a n (ll999h . where more details of the numeri- 
cal method are available. Typically, these techniques lead to an 
improvement of a factor ~ 2 - 3 in spatial resolution. 

The initial model consisted of an isothermal C-O white 
dwarf of 1.38 M0 in hydrostatic equilibrium, with a central den- 
sity pc = 1.8 X 10'' g cm"-*. We ran two simulations of the 
delayed detonation scenario, whose main features are given in 
Table [T] In model DDT3DA the first stages of the explosion 
were followed using a one-dimensional Lagrangian code until 
the mass burnt was 0.09 M0 and the central density had de- 
clined to pc = 1-4 X 10^ g cm""*. At that point the structure was 
mapped to a 3D distribution of 250 000 particles, and the radial 
velocity field at the flame location was perturbed by a sinusoid 
of wavenumber / = 20 and amplitude 200 km s"' , as much as the 
velocity of matter behind the flame at that time. The ensuing evo- 
lution was followed with the three-dimensional SPH code, with 
a baseline velocity of Vflame = 100 km s"'. The initial geometry 
of the flame in model DDT3DB was different, consisting in a set 
of 30 sparks scattered around the center of the white dwarf. The 
three-dimensio nal calculation of this model s tarted exactly as in 
model B30U of lGarcia-Senz & Bravol (l2005h . to which we refer 
the reader for more details. 

The switch from the subsonic combustion front to the super- 
sonic regime also was handled in a different way in both models. 
In DDT3DA the detonation was induced simultaneously in many 
points on the flame surface, when its average density declined to 
~ 2x 10^ g cm"^. By contrast, in model DDT3DB the detonation 
was artificially triggered in a single region of the flame, although 
at a density similar to that in DDT3DA. 

The last three columns of Table [T] give details of the per- 
formance of the code with respect to the conservation of total 
energy, momentum and angular momentum. The conservation 
of the two last is excellent, whereas that of energy is acceptable 
as the error stays around 1% of the final kinetic energy. Also 
supporting the ability of the code to handle the subsonic defla- 
gration phase is the spherical symmetry retained by the pressure 
field during the huge expansion of the white dwarf in this phase, 
as shown in Fig. [1] Such a smooth pressure distribution has to 
be compared with the convoluted temperature field depicted in 
the last snapshot of Fig. |2l The limitations of our methodology 
to carry out simulations of Type la sup ernovae were widely dis- 
cussed in lGarcia-Senz & Bravol (l2005h . 

3. Results of the three-dimensional simulations 

3.1. Model DDT3DA: The deflagration phase 

This phase begins when the first sparks ignite in the central re- 
gions of the white dwarf and ends when the first detonation 
wave appears. Both modes of burning propagation can coexist 
for a while during the detonation phase although the efliciency 
of the deflagration is far below that of the supersonic detonation. 
As mentioned before, in our computations the deflagration al- 




Fig. 1. Contour map of pressure in a cut across the star at the 
end of the deflagration phase of model DDT3DA (t = 1.55 s). 
The labels of the isobars give the logarithm of pressure in erg 
cm"^. The axis labels are in units of 1 000 km. The preservation 
of spherical symmetry is a consequence of the subsonic nature 
of the deflagration front. 



gor ithm was switch e d off o nce a DDT was artificially induced. 
Ro pke & Niemeveil (|200 7|) have discussed the impact of sub- 
sonic burning after a DDT. They found that it is relevant only 
when the mass burnt during the deflagration phase is quite low 
as is the case, for instance, if the runaway starts in a few hot 
bubbles. 

The geometry of the initial sparks which light the core 
of the white dwarf is not known. However, in the three- 
dimensional realm it is to be expected that the ignition starts 
asynchronously in a number of small regions scattered around 
the core (jOarcia-Senz & WooslevI 119951 : IWunsch & Wooslevi 
llOOl iKuhlen et alJl2006h . For a large enough number of ignit- 
ing points located close to the center, bubbles of incinerated mat- 
ter would merge before they have had the opportunity to detach 
from the burning pack. Were this the case, the hydrodynamical 
evolution would not be very different from that due to the mas- 
sive ignition of the central volume of the core. This is the sce- 
nario we examine in model DDT3DA, in which the explosion 
was triggered by incinerating a small amount of mass in the cen- 
tral part of the core of an isothermal Chandrasekhar-mass white 
dwarf. 

The evolution of the star during the deflagration phase can be 
seen in Figs.|2]and[3] In Fig.|2]there are several snapshots of the 
evolution from the beginning of the SPH calculation, which we 
define as f = s, when the flame surface was completely spher- 
ical, up to f = 1.55 s, when the flame was strongly distorted 
by the Rayleigh-Taylor instability and the average density at the 
flame was < poame >- 2 x 10^ g cm"-*. A complementary view 
is given in Fig. [3] where we plot the angle-averaged radial pro- 
files of density, temperature, nuclear energy generation rate, and 
incinerated mass fraction during the deflagration phase. During 
the first half a second the flame front remained nearly spherical 
in spite of the sample of perturbations of the velocity field that 
were seeded at f = s. Afterwards, the fingers characteristic of 
the Rayleigh-Taylor instability developed very quickly (second 
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Table 1. Main features of the models. 



Model 


Ignition 




M (5*Ni) 


flast 


/ 




^lf 


AL'' 






erg) 


(Mo) 


(s) 


(Time steps) 


(10^' erg) 


(lO'Mocm s-i) 


(10" M© cm^ s-') 


DDT3DA 


Central volume 


0.75 


0.54 


7.9 


7400 


0.01 


4 X 10-** 


2 X 10-** 


DDT3DB 


30 bubbles 


0.81 


0.63 


463. 


21200 


-0.008 


1 X 10-** 


10-** 



" Final energy (gravitational+internal+kinetic) of the ejected mass at the last computed time, rias,. 
Accumulated error in the total energy at the last computed model, after / time steps. 

Total momentum at the last computed model, normalized taking a characteristic velocity of the ejecta of 10' cm s"' . 

Total angular momentum at the last computed model, normalized taking a characteristic velocity of the ejecta of 10' cm s"' , and a characteristic 
radius of 10* cm. 



and third snapshots in Fig.|2]l, which resulted in hot material ris- 
ing and cool unburnt carbon and oxygen intruding towards the 
center (fourth panel of Fig. O. As the flame surface increased 
due to the hydrodynamic instabilities, the rate of released nu- 
clear energy also did. The effective burning velocity kept pace 
with the evolution of the surface, as can be seen in Fig. ID it re- 
mained around 0.02cs until t - 0.6 s and rose monotonically up 
to 0.14cs at f = 1.4 s. This large increase in the effective burning 
rate when t> 0.6 s is as well apparent in the third panel of Fig. [3] 
which depicts the nuclear energy generation rate. At f = 0.86 s 
(second snapshot in Fig.|2]i the burning front clearly displays the 
mushroom-like shape characteristic of the non-linear Rayleigh- 
Taylor regime, encompassing a large fraction of the white dwarf 
mass. As a consequence the rate of nuclear energy generation 
went up steeply until a maximum value of 1 .3 x 10''' erg s"' was 
achieved at t=l.l s. During the last second of deflagration prop- 
agation the average temperature and nuclear energy generation 
rate were uniform throughout the white dwarf. 

When the elapsed time was ~ 1 .55 s the average density of 
the flame became low enough to allow intermediate-mass ele- 
ments to be synthesized. At this time the central density had 
declined to pc - 5.4 x 10^ g cm"^ and the expansion was fast 
enough that the nuclear reactions began freezing at the outermost 
regions of the deflagration front. The main features of the model 
at the end of this phase are summarized in Table |2] The amount 
of nickel synthesized, 0.27 M0, was too low to give a suitable 
explosion with the maximum luminosity of a typical SNla. In 
addition, there remained too much carbon and oxygen. 

The scaling properties of the flame surface during the de- 
flagration phase were monitored, using the method described in 
[Garcia-Senz et al. ( 1998), and the results are presented in Fig.|5] 
(filled dots). This kind of calculation, which gives the so-called 
correlation dimension, provides a purely geometric description 
of the flame structure determined by the spatial distribution of 
particles belonging to the flame. An independent way of obtain- 
ing the fractal dimension is to estimate the total flame surface 
within the minimum lengthscale resolved by the code, h, and the 
integral scale, using the ratio of the mass consumption rate, M, 
and the baseline flame velocity (solid line in Fig. [Sj: 



D^2 + 



log [mi A: 



(1) 



log (rlh) 

where the integral scale has been approximated by the aver- 
age radius of the flame, r. As it can be seen in Fig. |5] there 
is excellent agreement between both methods of computation 
of the fractal dimension, implying that the above equation can 
be safely used in analytical and one-dimensional calculations to 
explore the physics of thermonuclear supernovae. After an ini- 
tial period where the fractal dimension remained close to 2.0 it 
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Fig. 3. Angle-averaged profiles of density, temperature, nuclear 
energy generation, and incinerated mass fraction during the de- 
flagration phase. The lines represent the profiles at the same 
times as Fig.|2l t - 0.56 s (soHd), 0.86 s(dotted), 1.20 s (short 
dashed), and 1.55 s (dot-long dashed). 



rapidly began to increase after t - 0.6 s, at which time the av- 
erage density of the flame had dropped below 5 x 10'^ g cm"^ . 
Approximately 1.1s after the beginning of the three-dimensional 
calculation the fractal dimension had reached a value D ^ 2.6, 
similar to that associated with the RT instability in the non linear 
regime, Drt = 2.5. That high value of the fractal dimension indi- 
cates that the Rayleigh-Taylor instability rather than turbulenc^ 
was the dominant mechanism wrinkling the flame at those scale- 
lengths directly resolved by the hydrocode. This was expected as 
our resolution is insufficient to resolve turbulence, i.e. the min- 
imum smoothing length attained in the SPH calculation during 
deflagration was 20 km while the integral length of the turbu- 



- In the d i scussion we are as suming Kolmogorov turbulence 
jZingale et alj 1 20051 ; ISchmidt et alj ]_2006), characterized by a frac- 
tal dimension Dk = 2.33, although the turbulence might also 
be of the Bolgiano-Obuko v type, characterized by £)bo = 2.6 
jNiemever & Kersteidl2000l) 
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Fig. 2. Snapshots of the temperature map in a cut across model DDT3DA at times t - 0.56, 0.86, 1.20, and 1.55 s, all of them 
belonging to the deflagration phase. The lengthscale is given in units of 1 000 km s while the temperature scale is given in the 
side colorbar in units of 10'' K. 



Table 2. Main features of the models at the time of DDT. 



Model 


DDT 


^DDT 

(s) 


^DDT 

(10' g 


(10' gcm-3) 


^tot 

(10=' erg) 


mDDT 

-'"nse 
(Mo) 


M(=^Ni)'^^T 
(Mo) 


M(C - 0)"^" 
(Mo) 


DDT3DA 
DDT3DB 


Multipoint 
Single point 


1.55 
0.81 


5.4 
22.2 


2.0 
3.7 


0.40 
0.39 


0.50 
0.58 


0.27 
0.43 


0.65 
0.64 



lent cascade is ~ 10 km dLisewski et al.ll2000l) . Afterwards, the 
fractal dimension saturated at D ^ 2.6 - 2.7. 

3.2. Model DDT3DA: The detonation phase 

As has been discussed in Section 2, explosions driven purely by 
a deflagration usually run into trouble when compared to the ob- 
servational properties of normal SNIa events for various reasons. 
One of them is the difficulty of reaching explosions with energies 



of roughly 10^' ergs. Another serious objection is that current 
three-dimensional calculations always predict a strong mixing 
of elements in the ejecta. A possible way to circumvent these 
drawbacks is to assume that at some time the subsonic deflagra- 
tion turns into a supersonic detonation at one or several points 
of the flame surface. Although several mechanisms have been 
suggested to account for such a transition, a definitive answer to 
this question is still missing. Therefore, the temporal and spatial 
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Fig. 7. Snapshots of the temperature map in a cut across model DDT3DA at times f = 1.58, 1.61,1 .66, and 1.75 s. The deflagration- 
to-detonation transition was induced just at 1.55 s. The scales of length and temperature are the same as in Fig.|2] 



coordinates at which the DDT takes place should be regarded as 
free parameters of the model (Livne 1999). 

Many one-dimensional simulations carried out in the 
nineties suggested that the agreement with observations is fairly 
good when the transition takes place at densities around 2 - 4 x 
10^ g cm~^. This is about the average density of the flame when 
its fractal dimension stabilizes around 2.6, and approaches the 
density at which combustion is thought to leave the flamelet 
regime and enter into the distributed regime. Let us assume that 
the flame propagates with the laminar velocity at the Gibson 
scale and that the flame surface behaves as a fractal of dimen- 
sion 2.33 (a value determined by the turbulent energy spec- 
trum) between the Gibson scale and the minimum scale resolved 
by the code, h. Then, the effective flame velocity can be ob- 
tained as a function of the fractal dimension, D, that charac- 



terizes the flame between h and the integral scale. The result 
of this calculation is that the effective flame velocity exceeds 
the maxim um velocity of a stable Chapman- Jouget deflagration 
(lKhokhlov.1988) if D > 2.5 atp $ 2 x 10^ g cm ^ Thus, even 
though the laminar flame velocity always remains much lower 
than the sound velocity, the effective flame velocity can be high 
enough to compress the material ahead, switching on a detona- 
tion ( IWooslev & WeaveiHl994b . 

As a practical criterion to select the particles prone to deto- 
nate we flagged those placed at regions where the local fractal 
dimension was higher than D = 2.5 when the average density 
of the flame was < poame >- 2 x 10^ g cm"-'. As these high-D 
regions trace the geometrical complexity of the flame it is not 
unreasonable to assume that they can harbor conditions suitable 
for detonation initiation. In our calculation the above criterion 
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Fig. 4. Effective velocity of the deflagration front in units of the Fig. 6. Flame structure in a cut of the white dwarf at the end of 
local sound speed as a function of time. the deflagration phase. The image shows the location of hot fuel 

whose temperature Ues in the range 1 - 3 x 10"^ K. 




duced through their instantaneous incineration. Note that we do 
not mean that reaching a given fractal dimension of the flame is 
a necessary or sufficient condition for transition to detonation, 
but for our purposes it provided a convenient way to select the 
time and location of the DDT. 



Fig. 5. Evolution of the fractal dimension of the flame. Dots 
show the correlation dimension, while the solid line represents 
an alternative estimate of the fractal dimension derived from the 
effective rate of mass burning. 



was satisfied for the first time at f = 1 .55 s, marking the end of 
the deflagration phase. A two-dimensional image of the flame 
structure at this time is shown in Fig. |6] where its geometric 
complexity can be seen. The regions displaying the highest D 
spanned the whole radial extent of the white dwarf, except for 
a narrow central volume and the outermost 0.2 M0. Once the 
high-D particles were selected, a detonation was artificially in- 
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Table 3. Nucleosynthesis products, in solar masses. 



0.4 - 




Fig. 8. Evolution of angle-averag ed profiles of the i^c + i^O 
mass fraction during the detonation phase of model DDT3DA. 
The profile at DDT is shown together with those belonging to 
the same times shown in Fig.|7] 



As we did with the deflagration phase, we describe here the 
evolution during the detonation phase by means of two figures. 
Figs. |2] and [8] the first one showing several snapshots of the tem- 
perature field through a cut of the white dwarf, and the second 
one giving the evolution of the angle-averaged radial profiles of 
flie fuel mass fraction. After a transient induction period, deto- 
nations settle in a steady state consisting of a shock plus a tiny 
reaction zone, followed by a much larger relaxation zone. Our 
code is only able to effectively resolve the last one. Thus, we 
cannot ascertain the result of the interaction of the detonation 
front with obstacles, i.e. inert volumes. However, matter inciner- 
ated to the NSE should not be described as inert, because at the 
densities we are considering it is able to react to the passage of 
a shock, that compres s and heats matter, and re adjust its com- 
position accordingly jMaier & Niemeve^ l2006h . Hence, when 
the precursor shock of a detonation transits through NSE mat- 
ter, it rises its temperature favouring lighter, less bound, species, 
which in turn depletes energy from the shock front. In our cal- 
culations, the detonation did not survive the crossing of volumes 
composed of nuclear ashes. The complex geometry of the flame 
at the time of DDT favored the formation of isolated pockets of 
fuel that could not be caught by the detonation front, so that they 
remained unburned (Fig.|9j. 

Because the detonation was initiated in a more or less uni- 
form way throughout the white dwarf, there followed a period 
of rapid and generalized combustion leading to a temperature 
distribution much more homogeneus than that of the preced- 
ing deflagration phase (Fig. |2l). The inwards detonation wave 
synthesized an additional mass of ^^Ni, whereas the outgoing 
detonation, travelling through lower density layers, synthesized 
intermediate-mass elements (Fig. |9]l. In this way, the DDT al- 
lowed the partial compensation of the shortcomings of the pre- 
vious deflagration stage, giving an improved explosion model. 
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Fig. 9. Map of the final chemical composition on a cut across 
the ejected matter: C-O is shown in black, intermediate-mass el- 
ements in grey, and Fe-Ni in white. The incinerated matter per- 
vades the whole ejecta, with intermediate-mass elements domi- 
nating an annular region at intermediate radius. Unburned C-O 
are found in isolated pockets, mainly in the outer regions of the 
ejecta. The origin of the large C-O pockets seen in this figure 
can be traced back to features of the flame front identifiable in 
Figs.Eliandll] 



The main features of the explosion at the last computed 
time are given in Table [1] while the nucleosynthesis obtained 
by post-processing the thermal history obtained in the hydro- 
dynamical simulation is presented in Table |3] The final kinetic 
energy is still on the low side but no longer very low whereas the 
amount of ^^Ni (0.54 Mq) is enough to power the light curve. 
Nevertheless, there still remains 0.39 Mq of unburnt carbon and 
oxygen, which should be compared with the 0.65 Mq left at 
the end of the deflagration phase. The reason why there still re- 
mains this large mass of carbon and oxygen is due to both the 
shielding effect of the ashes produced during the deflagration 
phase and the degradation of the code resolution when the det- 
onation propagated through low-density regions. However, even 
if a much better resolution were achieved we would not expect 
major changes in the total energy of the ejecta or in the mass of 
^''Ni synthesized. 
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Fig. 10. Final distribution of elements in model DDT3DA as a 
function of velocity, after radioactive decay. 



The total amount of Fe (after radioactive decay) and 
intermediate-mass elements ejected in the explosion are 
0.66 M0 and 0.11 Mq respectively, which is of the same order 
as the abundances deduced from the spectra of normal Type la 
supernovae and compatible with the demands of current models 
of galactic chemical evolution. Another positive feature of the 
model is that the ashes are not concentrated in large clumps as is 
usually obtained in pure deflagration models. In addition, there 
remains little unburnt carbon and oxygen close to the center. 
Unfortunately, model DDT3DA does not completely remove all 
nucleosynthetic deficiencies affecting current three-dimensional 
models of SNIa. There are too many Fe-group elements close 
to the surface of the white dwarf, especially too much radioac- 
tive ^^Ni (Fig. |9]l. The distribution of species in velocity space 
(Fig. [Toll shows that the ejecta is not chemically stratified, as de- 
manded by recents observations of the light curve in the near and 
middle infrared (Hoflich et al. 2002) as well as by the properties 
of the X-ray em ission of putative Type la supernova remnants 
jBadenes et al.ll200 6). The Fe that makes up most of the ejecta 
(Fig. [TOl see also Fig.fTTTi is mostly produced during the explo- 
sion in the form of radioactive ^^Ni, whereas the curves labelled 
as 'Ni' trace the distribution of stable iron-group species. Such 
stable nuclei are nearly absent from the low velocity layers of 
the ejecta, which doe s not match the requirements from infrared 
observatioris of SN Ia dMotohara et al.l2006l:lGerardv et alj2007t 
iFesen et alj|2007l) . 

3.3. Sensitivity studies 

We have explored the sensitivity of our results to the parameters 
of the initial model and to the deflagration-to-detonation transi- 
tion criterion. 

To start with, we changed the wavenumber of the perturba- 
tion imposed to the velocity field at the flame location at the 
beginning of the three-dimensional calculation. Decreasing the 



wavenumber from / = 20 to / = 7 has virtually no eff^ect on 
the results of the simulation, as the global figures change by less 
than 2-3%. Next, we checked the sensitivity to the amplitude 
of the perturbation. Although this parameter is more influential 
than the wavenumber, its imprint is still quite modest. Reducing 
the amplitude from 200 to 100 km s"' translates into a less ener- 
getic explosion (0.16 foes less than in DDT3DA) and produces 
0.11 M0 less of ^''Ni. Of course, this implies that the remaining 
mass of unburnt carbon and oxygen is larger than in DDT3DA. 

To explore the sensitivity of the results to the DDT criterion 
we took the same configuration at the end of the deflagration 
phase as in DDT3DA, but used alternative prescriptions to select 
the location of the induced detonation: 

- Core detonation. All flame particles whose radius was lower 
than 1.3x10^ cm were instantaneously incinerated. The den- 
sity in these regions was higher than ~ 4 x 10^ g cm"^. 

- Mid-altitude detonation. In this variant all flame particles 
within radii 1.7 x 10** cm and 2.0 x 10** cm were instanta- 
neously incinerated. Their mean density was approximately 
~ 2 X 10^ g cm-\ 

- Atmospheric detonation. All flame particles within 2.4 x 
10** cm and 2.8 x 10** cm were instantaneously inciner- 
ated. Their densities ranged from ~ 7 x 10^ g cm"^ to 
~ 10^ g cm'\ 

The diversity of supernova explosions obtained with these 
prescriptions is not as rich as that usually achieved in one- 
dimensional delayed detonations. The kinetic energy varies at 
most in 0.15 foes, whereas the range of ""^Ni masses spans from 
0.34 to 0.43 M0. By contrast, one-dimensional models with 
transition densities in the same range as above produce explo- 
sions with kinetic energies that change by 0.34 foes and ^^Ni 
masses rangmg from 0.38 to 0.97 M© (Bad enes etaIll2003L 
2005). Other properties of these three-dimensional explosions 
are similar to what was found in DDT3DA: the remaining mass 
of unburned carbon and oxygen is in all cases too large, and they 
lack the desired level of chemical stratification. 

3.4. Model DDT3DB 

The evolution of model DDT3DB during the first 0.81 s belong- 
in g to the deflagration phase wa s identical to that of model B30U 
in iGarcfa-Senz & Bravol (l2005h . At fliat time a DDT was artifi- 
cially induced in the vicinity of the flame. The properties of the 
model at DDT are given in Table |2] In comparison with the con- 
figuration of DDT3DA at the moment of transition to detona- 
tion, the structure of DDT3DB was more centrally condensed, 
the central density being a factor ~ 4 higher and the average 
flame density nearly 3.7 x 10^ g cm"^. The total energy and in- 
cinerated mass are, however, quite similar in both models, as is 
the amount of unburned fuel. The major difference was in the 
mass of ^^Ni synthesized during the deflagration phase, due to 
the longest time spent by DDT3DA at high densities, which en- 
abled incinerated matter to experience more electron captures. 
Consequently, the deflagration phase of model DDT3DB was 
able to synthesize more ^^Ni, but less neutronized elements. This 
difference was still visible in the final results, after the detonation 
phase, although somewhat attenuated (Table [TJ. 

The temporal evolution of the detonation phase in model 
DDT3DB is shown in Fig. [TTI As can be seen in Fig. 5 of 
iOarcfa-Senz & Bravol (l2005l central figure in the right column) 

^ Note that, at f = 0, the wavelength of the perturbation was ~ 5 times 
the resolution of the code for our reference model with / = 20 
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from first physical principles, this is an exceedingly difficult task 
in view of the many different length and temporal scales in- 
volved and of the huge diversity of initial conditions that can be 



Velocity (1,000 km/s) 



Fig. 12. Final distribution of elements in model DDT3DB as a 
function of velocity, after radioactive decay. 



at the time of DDT the hot nuclear ashes formed a broken ring 
that enclosed a central volume rich in carbon and oxygen. The 
ensuing detonation preferentially propagated inwards processing 
most of the inner fuel into Fe-peak nuclei without having to deal 
with major obstacles composed of almost inert matter. However, 
a detonation wave was not successfully launched outwards and, 
as a result, the already vigorous expansion of the white dwarf 
prevented the outer layers from being processed by the nuclear 
bomb. 

Model DDT3DB shares the same problems as DDT3DA: 
low final kinetic energy, large mass of unburned carbon and oxy- 
gen, and lack of chemical stratification (Fig. \T2\ . A significa- 
tive difference between both delayed detonation models is the 
amount of neutronized Fe-peak nuclei ejected in the explosion. 
The mass of slightly neutronized nuclei such as ^"^Fe and ^^Ni 
was of the same order in both models because their production 
is determined by the initial presence of ^"Ne ( 1 % in mass) mixed 
with '^C and '^O, rather than by electron captures on NSE mat- 
ter. The final ratio of the masses of ^^Fe and ^**Ni with respect 
to ^^Fe is about 1.5 and 2 times solar respectively. However, 
the production of more neutronized Fe-peak nuclei in model 
DDT3DB is more scarce than in DDT3DA due to the rapid mi- 
gration of the bubbles to low density regions shortly after the be- 
ginning of the deflagration phase. Had our initial models started 
from even higher central densities this effect would had been ac- 
centuated. 



4. Conclusions 

The multidimensional delayed detonation model still remains 
an alternative way to account for the observational properties 
of Type la supernovae, although there still remain many pieces 
of the puzzle to be settled. Among them is the mechanism of 
transition from deflagration to detonation, especially its viabil- 
ity in an exploding white dwarf. Although the most satisfac- 
tory approach to the problem would be to reveal the mechanism 



(iKhokhlov 1991 


i Khokhlov et al.' '19971; 'Niemeyer & Wooskvl 


19971 Niemevei 


1999; Lisewski et al. 2000; Ropke et al. 200l 


Zineale & Dursi 


2007|) no sound mechanism for the transition 



has yet been found. 

Given the lack of a consistent physical picture of the DDT 
mechanism the opposite path is usually taken, trying to con- 
strain the conditions of DDT by comparing the resulting explo- 
sion with the observational data of SNla. However, this approach 
is encumbered by the vast set of DDT conditions that need to be 
tested if a three-dimensional model has to be built. In this pa- 
per, we have tried to contribute to the still scarce set of three- 
dimensional simulations of delayed detonation models with two 
simulations in which the transition takes place at densities in the 
same range as those favored by succesful one-dimensional de- 
layed detonation models of SNIa, i.e. pti ~ 1 - 4 X 10^ g cm~^. 

Our delayed detonation models improve on the pure defla- 
gration ones, but there still remain severe deficiencies compared 
to observations of SNIa, the two most important being the lack 
of chemical stratification and the huge mass of unburned car- 
bon and oxygen ejected in the supernova explosion. While the 
amount of ejected carbon and oxygen might be affected by reso- 
lution issues, it is unclear if it can be lowered as much as needed 
to fit the observations. On the other hand, the lack of chemical 
stratification may not depend as much on the details of the nu- 
merical implementation. To our knowledge, a chemically strati- 
fied three-dimensional d elayed detonation has only been found 
in a single work so far (iGamezo et al.l l2005l) . but they had to 
assume a singularly high density of DDT. The presence of high- 
velocity nickel above the silicon layer could also be at odds with 
the observations of SNIa, but this cannot be ascertained in the 
absence of a three-dimensional spectral calculation. 

Another interesting result is the relative insensitivity of the 
global properties of the explosion with respect to variations in 
the initial model (geometry of the runaway volume: centered 
vs bubbly) and location of the DDT (at a fixed transition time). 
Differences in kinetic energy of 0.15 foes and ejected masses of 
^''Ni of 0.09 Mq as we obtained are in clear contrast to the rich 
diversity of explosions obtained with one-dimensional and two- 
dimensional (Arnett & Livne 1994) delayed detonation models 
in which the DDT spanned similar density ranges. 

Model DDT3DB can be compare d with the simulations per- 
formed by 'Ro pke & Niemeveij(l2007l) . In this last work three de- 
layed detonation models were presented in which the initial run- 
away took place in different numbers of hot bubbles, whereas the 
densities of transition to detonation ranged from 1 .3 to 2.4x 10^ g 
cm"^, compared with the transition density of model DDT3DB, 
Ptr = 3.7 X 10^ g cm"^ (Table|2]i. In spite of having used diverse 
numerical techniques and slightly different physics in both cal- 
culations, the final kinetic energy of our model DDT3DB (0.81 
foes) fits quite well in the trend of the asymptotic kinetic ener- 
gies given by Ropke & Niemeyer (1.5 to 1.0 foes). The same 
can be said about the ejected mass of unburned carbon and oxy- 
gen, which goes from 0.04 to 0.22 Mq in their models, to be 
compared with 0.45 M© of C-O ejected in DDT3DB. The con- 
vergence of results of these different simulations supports the 
models reported here. 
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